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It is often impractical to instrument the external surface of high-speed vehicles due to the 
aerothermodynamic heating. Temperatures can instead be measured internal to the structure 
using embedded thermocouples, and direct and inverse methods can then be used to estimate 
temperature and heat flux on the external surface. Two thermocouples embedded at different 
depths are required to solve direct and inverse problems, and filtering schemes are used to 
reduce noise in the measured data. Accuracy in the estimated surface temperature and heat 
flux is dependent on several factors. Factors include the thermocouple location through the 
thickness of a material, the sensitivity of the surface solution to the error in the specified 
location of the embedded thermocouples, and the sensitivity to the error in thermocouple data. 
The effect of these factors on solution accuracy is studied using the methodology discussed in 
the work of Pizzo, et. al.! A numerical study is performed to determine if there is an optimal 
depth at which to embed one thermocouple through the thickness of a material assuming that 
a second thermocouple is installed on the back face. Solution accuracy will be discussed for a 
range of embedded thermocouple depths. Moreover, the sensitivity of the surface solution to 
(a) the error in the specified location of the embedded thermocouple and to (b) the error in 
the thermocouple data are quantified using numerical simulation, and the results are 


discussed. 
Nomenclature 
1-D = One-dimensional 
Cp = Specific heat capacity, J/(kg-K) 
DHCP, IHCP = Direct (respectively, inverse) heat conduction problem 
€,(t),er(t) = Error in heat flux (respectively, temperature), % 
k = Thermal conductivity, W/(m-K) or W/(cm-K) 
L = Thickness of material specimen 
m = Superscript denoting the time step count 
n = Subscript denoting the nodal count 
Np, N, = Number of nodes for the DHCP (respectively, IHCP) 
q = Heat flux, W/m? or W/cm? 
T = Temperature, K 
t = Time, s 
Xn = Spatial location, m or cm 
a = Thermal diffusivity = k/pcp, m’/s 
At = Time step size, s 
Ax = Cell size, mor cm 
€ = Symbol denoting the phrase is an element of 
p = Mass density, kg/m? or kg/cm? 
* Ph.D. Candidate, Pathways Student Trainee 
> ATAA Associate Fellow 
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I. Introduction 

hile traveling at hypersonic speeds, flight vehicles encounter high heating resulting in high surface temperature. 

There is an interest in knowing these surface values particularly in the fields of materials, structures, and 
aerothermodynamics research. Unfortunately, it is usually either impractical or impossible to instrument the external 
surface of high-speed vehicles due to the aerothermodynamic heating. Temperature histories are instead measured 
internally to the surface by embedding thermocouples or thermocouple plugs through the thickness of the vehicle skin 
or thermal protection system. The measurements can then be used to solve an inverse heat conduction problem (IHCP) 
to estimate surface temperature and heat flux. The inverse problem is ill-posed. Therefore, the solution is extremely 
sensitive to measurement errors, i.e., small errors in the data yield large errors in the solution*. 

The study of inverse heat conduction started in the late 1950s with work published by Shumakov’ and Stolz® in 
1957 and 1960, respectively. In the years since, numerous methods for solving both linear (constant thermophysical 
properties) and nonlinear (temperature-dependent thermophysical properties) one-dimensional (1-D) transient inverse 
problems have emerged from researchers in the field such as Beck*?!°!!, Murio®, Weber’, and Carasso°, among others. 

In 2016, Pizzo, et. al.,! estimated surface temperature and heat flux histories of high-temperature carbon/carbon 
materials using a 1-D inverse method (defined within the body of paper) with temperature-dependent thermophysical 
properties (nonlinear problem) assuming negligible in-plane temperature gradients. The study was motivated by 
Frankel, et. al.,*"'?-'3. In his work, Frankel solved the 1-D transient heat conduction problem using noisy data measured 
at internal depths x, and x2 where x; < x2. He reduced noise in the data with a global low-pass Gaussian filter by 
defining a cutoff frequency using discrete Fourier transforms of the data. Frankel then approximated the 
heating/cooling rate using finite differences and numerically integrated the integral relationship between temperature 
and heat flux to obtain an analytical heat flux solution at depth x,. The surface temperature and heat flux were 
approximated from the measured temperature and calculated heat flux at depth x2 using finite differences and Taylor 
series reconstruction. This method yielded accurate results. 

The methodology in Pizzo, et. al' used Frankel’s procedures as a guide, but made two slight modifications. First, 
local windowed-sinc filters were used to reduce noise in the measured data. Second, a technique developed by Carasso° 
was used to solve the inverse problem. In this technique, the 1-D heat conduction equation is first written as a first- 
order system with temperature and heat flux as the unknowns, and is then discretized using forward differencing in 
space and central differencing in time. The methodology in [1] was validated using data provided by thermal vacuum 
chamber radiant heating tests (c.f. Blosser'*). Pizzo demonstrated accuracy in estimating temperature histories using 
comparisons with test results (estimated solutions had error in the range of -5% to 9%). The estimated heat flux 
histories however had significantly larger error than temperature (estimated solutions had error in the range of -30% 
to 10%). The larger errors were expected, because for the methodology used, the rate of convergence with heat flux 
with mesh size is one order less accurate than temperature. 

The purpose of this paper is to expand upon the work in [1] by characterizing the selected method for inverse heat 
conduction. For inverse methods, because small errors in the data yield large errors in the solution, accuracy in the 
estimated surface temperature and heat flux is dependent on several factors. Factors include the embedded 
thermocouple depth through the thickness of a material, the sensitivity of the surface solution to the error in the 
specified location of the embedded thermocouple, and the sensitivity of the solution to the error in thermocouple data. 
The effect of these factors on solution accuracy is presented in this paper. 

The inverse methodology is summarized in Section I, and the characterization of the method is discussed in 
Section III. Concluding remarks are discussed in Section IV. 


II. Methodology Summary 


The transient temperature distribution in a material is governed by the heat conduction equation. Assuming that 
the in-plane temperature gradients are negligible, the heat conduction equation reduces to its 1-D form: 


ar a aT 
Pep (T) Fe = ay (cr) ak . 


where p, C,, and k are the density, specific heat, and thermal conductivity of the material, respectively. The initial 
temperature distribution through the thickness of the material is specified at time t = 0: 


T(x, 0) = T;(x). (2) 
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Boundary conditions at the outer surfaces (x = 0, L) are required to solve Eq. (1). For the classical initial-boundary 
value problem, referred to in this work as the direct heat conduction problem (DHCP), the boundary conditions are in 
the form of specified temperature or heat flux*. In situations where the boundary values are unknown, an inverse**!” 
problem must be solved using known internal 
temperatures at two depths, denoted by x, and x2. The (L, t) 
source of the internal temperatures is typically data od is cio a TERR ee enn x =L 
from a set of embedded thermocouples as illustrated ace, | = 
in Figure 1. 

In [1], the surface temperature and heat flUX  ssrretssseteseeteeeessees 
histories at x = L are estimated from two internal 
temperature measurements, 7, and T>, using a four- 
step process. The four-step process includes filtering 


Heat conducting 
material 


noise from the measured data, solving for temperature pas 
histories in the region between depths x, and x2 
where 0 < x1 < x2 <L, solving for the heat flux = 0 


history at depth, x2, and solving for the temperature Figure 1. Illustration of the direct and inverse problems 
and heat flux histories in the region extending from and the embedded thermocouples, T: and T2. 
depth x =x, tox =L. 


Solving for Temperature Histories Between Depths x; and x, (Direct Problem) 
The region between thermocouples embedded at depths x, and x, is treated as the DHCP with initial and boundary 
values given by measured data, i.e. 


T(x,0) = Tp, XS xXx 

3a 
T(x,t)=T,  t>0 GB) 
T(%,t)=T,  t>0. 


Equation (1) is discretized by dividing the domain into a uniform mesh of Np total nodes with cell size of Ax = 
(x2 — x,)/(Np — 1). In the mesh, node n = 1 is at depth x = x, and node n = Np is at depth x = x,. Letting 7,” 
denote the temperature of node n at time t,, = mAt, the temperature of node n at time ty,4, = (m + I)At is obtained 
using a vertex-based finite volume method with Crank-Nicholson time marching: 


m+1 m m m 
Th ~ Th _ An-1)0n Annt+1 


+1 1 
At = 2Ax2 Git, Pie coal )+ 2Ax2 


(Tita — Ta” + Trea” — Tr?) (4) 


where 7 14.n (Anin+1) denotes the average of the temperature-dependent thermal diffusivity evaluated at temperatures 
Trt, and T,” (Tz” and T;,1). Since T]"** = T(x,,(m + 1)At) and Tyt** = T(x2,(m + 1)At) are known from the 
filtered data, Eq. (4) yields temperature histories at nodes n € [2, Np — 1], or equivalently, in the domain x € [x,, x2]. 


Solving for the Heat Flux History at Depth x, (Direct Problem) 

Let qn, denote the heat flux entering the domain at node Np at time t,, = mAt. The heat flux entering the domain 
at node Np at time t,,4, = (m+ 1)At is treated as an unknown and is obtained using a consistent approach based 
upon an energy balance for the half-call x € [xvp — Ax/2, Xn | combined with Crank-Nicholson time marching: 


Node Np: Qn’ = ki, 


TN Np One oo ‘ie pect (2) ae - “i is 


2x “Pnp\ 9 At 


where ky, and Con, denote the temperature-dependent thermal conductivity and specific heat, respectively, evaluated 
D 


at temperature Ty. The density is denoted by p and is constant. Solving Eq. (5) yields the heat flux history at node 
n = Np, or equivalently, at depth x = x. 
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Solving for Temperature Heat Flux Histories Between Depths x2 and x = L (Inverse Problem) 

The region extending from thermocouple depth x = x, to the external surface x = L is treated as the IHCP. The 
temperature and heat flux distributions are calculated in the domain x € [x2, L] using a mesh of N; nodes where node 
n = 1 is at depth x = x2 and node n = N, is at x = L. The method used to solve the inverse problem is the space 
marching scheme S6 as outlined by Carasso>. The scheme has a truncation error of O(At?) and O(Ax), and the scheme 
was found by Carasso to have the smallest amplification factor. Therefore, the S6 scheme was the most stable of the 
O(At?) schemes analyzed. The S6 scheme is derived by writing the 1-D heat conduction equation as a first-order 
system for the heat flux: 


aT aq aT 
ae ax’ ax P® ce) ot’ 


The temperature and heat flux at node 1 are known since x, is a shared boundary between the domain of the direct 
and inverse problems. Equation (6) is discretized using forward differencing in space and central differencing in time. 
For each node, n = 1,..., N; — 1: 


Ta =T,; and Gnet = dn (7) 


Tria = Tr + Ax(qr/ke) and qny1 = Qn + Axpcp™ (Tat? —T,"*)/24t, m=1,..,M—1 (8) 


Ta = oH Ba a Tee and Gra = 204 ay Aner (9) 


Equations (7) to (9) yield temperature and heat flux distributions in the domain x € [x,,L] marching from node 1 
to node N, using boundary values given by the DHCP solution at depth x9, i.e. 


T (x2,¢) = Tpucr, x, and q(X2,t) = Qpucp, x, (10) 


where Tpycp, x, denotes the filtered temperature at depth x, and qpycp, x, denotes the heat flux solution at depth x, 
calculated from Eq. (5). Starting with node 1, the solution is marched in time before preceding to the next node. 


III. Methodology Characterization 


The methodology in Pizzo, et. al,' is characterized by studying (a) the effects of embedded thermocouple depth 
through the thickness of a material on the accuracy of the estimated surface temperature and heat flux, (b) the 
sensitivity of the surface solution to the error in the specified location of the embedded thermocouples, and (c) the 
sensitivity to the error in the thermocouple data. For the characterization, it is assumed that one thermocouple is 
installed on the back face of a material and a second thermocouple is embedded through the thickness. The motivation 
for this work is to study the factors that affect solution accuracy and identify best practices when solving inverse 
problems using the methodology presented. 


Model Problem 
Each numerical study is performed using a model problem for which an analytical solution exists. The model 
assumes that one thermocouple is embedded through the thickness of a material and a second thermocouple is installed 
on the back face. Numerical solutions are compared to the exact analytical solution to assess accuracy of the method. 
Assuming constant thermophysical properties, the heat conduction equation (1) reduces to: 


eve ae (11) 


where a = k/pcp. In each study, Eq. (11) is solved analytically in 0 < x < L for t € [0,100] s, with initial and 
boundary conditions: 
T(x, 0) = 294K, O< x<L 


k (OT /0x)(0,t) = qo, t>0 (12) 


k (OT /0x)(L,t) = qz, t>0 
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and thermophysical properties: 


k = 0.0575 Wm-K),  p = 0.0015 kg/cm, cp = 1833.33 J(kg-K). (13) 


Boundary conditions of qo = 0 W/cm?, q, = 40 W/cm? (herein identified as the Constant Case) and qg = 0 W/cm?, 
q, = (40 — t/10) W/cm? (herein identified as the Linear Case) are considered in Eq. (12). The imposed heat flux 
boundary conditions at depth x, are graphed versus time in Figure 2 for both the (a) Constant Case and (b) Linear 
Case, where t € [0,100] s. The corresponding analytical temperature solutions at depth x,, denoted by T,, are graphed 
versus time in Figure 3 for both the (a) Constant Case and (b) Linear Case. 


50 % 50 2000 2000 
a 2 
5 40 2 | 1500 1500 
e = = = fe 1000 g 1000 
= an 8 
$ 20 = 20 ee wo 
wv 
= 10 = 10 500 500 
=| 
% 25 50 75 100 e % 25 50 75 100 % 50 100 % 50 100 
Time, s Time, s Time, s Time, s 
(a) Constant Case (b) Linear Case (a) Constant Case (b) Linear Case 
Figure 2. Illustration of the imposed heat flux q, Figure 3. Illustration of the analytical temperature 
graphed versus time. T, graphed versus time. 


For each of the numerical studies presented in this paper, solution convergence was achieved by refining the spatial 
domain of the direct problem. With each successive refinement, the cell size was divided in half and the unbiased 
variance was calculated between consecutive refinements using the temperature of the node at the midpoint in the 
spatial domain. Once the variance fell below 0.01 K, the solution was considered converged using the refinement of 
the smaller cell size. Moreover, time-steps of At = 0.1 s and At = 0.01 s were considered in each study. 

The estimated and exact temperature at x = L are used to compute the error in temperature: 


(14) 


T, 


estimated T(L, t) — T, 
et =100| —————— 


where the exact temperature T, is obtained analytically. The graphical representations of the analytical temperatures 
are shown in Figure 3. The estimated and imposed heat flux at x = L are used to compute the error in heat flux: 


(15) 


se T80 a q(L, t) - 1) 
qd 


QL 


where the imposed heat flux is given by q, = 40 W/cm’ in the Constant Case and q, = (40 — t/10) W/cm? in the 
Linear Case. The graphical representations of the imposed heat flux are shown in Figure 2. The absolute values were 
removed from the standard error calculation in Eqs. (14) and (15) to examine when the numerical solution either 
overestimated (e > 0) or underestimated (e < 0) the solution. 


Effects of Embedded Thermocouple Depth Through the Thickness of a Material 

A numerical study was conducted to determine if there is an optimal depth to embed one thermocouple through 
the thickness of a material. Fifteen different depths were considered for the location of the embedded thermocouple, 
namely depths X,,Xz,...X% 14, 5 as illustrated in Figure 4. Using the analytical temperature solutions of Eqs. (11) to 
(13) at depth x9 paired with the solutions at depths x1, X2,...X14,X15 aS inputs to the DHCP, heat flux histories are 
estimated at depths x1, X2, ...X14,%,5. The boundary values in Eq. (10) are then given by the analytical temperature 
and estimated heat flux at each depth, and are used to solve the IHCP to estimate the surface temperature and heat flux 
histories, T;, and q,, respectively. 
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k (OT /0x) (L,t) = — es Xf) | L 


X15 = 15L/16 
Se ee eee ee eee ee X44 = 7L/8 
2 ee eee eee fhe, X13 = 13L/16 

_ X42 = 3L/4 
Bo rn er eo ee er Eerie) eet er eet X14, = 11L/16 
22-------- +--+ +--+ ~~ - == fh =} --f --b--. X19 = SL/8 

Xy =9L/16 
22---------------------- --- f -- f --  --} --f --f --b--. Xg =L/2 
-----------------------}--} --4} --} --} --} --}--}--. x7 =7L/16 

X_ =3L/8 
aoe eee ff -- f -- 4} -- } -- f fb bf = x, = 51/16 
eoee eee ff dd hb hp bbe x, = 1/4 

x3 =3L/16 
--------4--4--4--4--4--4--4--4--} --} --} --L--}--. x, = 1/8 
ng ee nnn dn dn nnn ncn 4 = L/16 

k (OT /0x) (0,t) = qo Xo =0 


Figure 4. Illustration of the fifteen different depths considered for the embedded thermocouple. 


In both the Constant and Linear Case studies with time-step of At = 0.01 s, the ill-posedness of the inverse 
problem resulted in growing oscillations as the solution marched from each simulated thermocouple depth to x,. The 
oscillatory growth was greatest for depths x, to x, yielding meaningless results with embedded thermocouples nearest 
the back face. The oscillations are a direct result of the domain size in the inverse problem combined with a small 
time-step. With small time-steps, larger domain sizes for the inverse problem allow for greater oscillatory growth. On 
the other hand, in both the Constant and Linear Case studies with time-step of At = 0.1 s, the effect of the ill- 
posedness of the inverse problem was negligible for all embedded thermocouple depths x, to x,;. Temperature and 
heat flux estimations are therefore only given for the Constant and Linear Case studies with time-step of At = 0.1 s. 
Moreover, only the estimations from depths x7 to x,, are considered for assessment due to the growing oscillations 
observed for depths x, to x, with the smaller time-step of At = 0.01 s. 

The Constant Case temperature and heat flux estimates are shown in Figure 5(a) and Figure 5(b) for At = 0.1 s. 
Neglecting depths x, to xg, the smallest magnitude of error in temperature is obtained for depths Xg, x11, X12, X7, and 
X19 with error ranging from 0% to 4% between t € [0,10] s, and from 0% to 2% by t = 100 s. The largest magnitude 
of error in temperature is obtained for depths x13, x15, X9, and x 4 with error greater than 5% between t € [0,10] s, 
and from 2% to 3% by t = 100s. The smallest magnitude of error in heat flux is obtained for depths x12, x11, X19, X7, 
X13, and x,, with error ranging from —1% to 1% for t € [0,100] s. The largest magnitude of error in heat flux is 
obtained for depths x9, xg, and x,, with error greater than +2% for t € [0,100] s. 

The Linear Case temperature and heat flux estimates are shown in Figure 6(a) and Figure 6(b) for At = 0.1 s. 
Neglecting depths x, to x, the smallest magnitude of error in both temperature and heat flux is obtained for depths 
X15, X44, X13, X12, X41, and xq with error decreasing below —0.2% by t = 100 s for temperature, and decreasing 
below 0.05% by t = 100 s for heat flux. 

In the Linear Case, depths x13 to x;5 yield the smallest magnitude of error in both temperature and heat flux, 
followed closely by depths x19 to x,2. However, in the Constant Case, depths x3 to x1, yield the largest magnitude 
of error in both temperature and heat flux, whereas depths x19 to X,2 yield the smallest. Considering all of the 
numerical results, the studies demonstrate that when using the methodology presented in [1] the surface estimations 
for both temperature and heat flux are consistently the most accurate if one thermocouple is embedded between depths 
X49 and X12 (5L/8 and 3L/4) when the second thermocouple is installed on the back face. This finding is based entirely 
off of the model problem previously discussed, and is independent of time-step. Further assessments should be 
conducted to conclusively make this assertion, as alluded to in the concluding remarks in Section IV. 
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(b) Heat flux error at x = L 


Figure 5. Error in the Constant Case estimations using a time step of At = 0.1 s. 
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(b) Heat flux error at x = L 
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Figure 6. Error in the Linear Case estimations, comparing the numerical results to the analytical solution with 


considering time-step of At = 0.1 s. 
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Sensitivity of Surface Solution to Error in Specified Location of the Embedded Thermocouple 

A numerical study was conducted to assess the sensitivity of the surface solution to the error in the specified 
location of the embedded thermocouple. For the sensitivity study, the embedded thermocouple is assumed to be at 
depth 3L/4. The depth of 3L/4 was selected because the previous results demonstrate that the estimations for both 
temperature and heat flux are consistently the most accurate when the embedded thermocouple is between depths 
5L/8 and 3L/4. 

To simulate error in the specified location of the thermocouple, the stated methodology is solved with accurate 
temperature data located at inaccurate depths. The embedded depth is perturbed by a factor of 6. For 6 = 0.01, 0.02, 
0.03, 0.04, 0.05, 0.10, and 0.15, the DHCP is solved with analytical solutions of T(x,,t) and T (xz + 6L, t) and depth 
inputs of x; = 0 and x, = 3L/4. The solution yields qpycp, x, defined by Eq. (10). Boundary values T(x. + 6L,t) and 
Apucp, x, are then used to solve the IHCP to estimate the surface histories, T, and q,. 

Surface estimations are compared against the error in temperature and heat flux, e and e, respectively, obtained 
from the embedded thermocouple depth study. The known e; and e, from depth x,2, shown in Figure 5 and Figure 
6, are used as a baseline in which case the perturbation, 6, is equal to 0.00. Since the model problem is a linear problem, 
the difference in error in the estimated temperature and heat flux is expected to be proportional to the prescribed error 
in the simulated baseline data. 

The Constant Case temperature and heat flux estimates are shown in Figure 7(a) and Figure 7(b) for At = 0.1 s. 
Additionally, the Linear Case temperature and heat flux estimates are shown in Figure 8(a) and Figure 8(b) for At = 
0.1 s. Results demonstrate that for both the Constant and Linear Case, as 6 increases (decreases) towards +0.15 
(—0.15), the error in surface temperature and heat flux increases (decreases). Hence, the difference in error in the 
estimated temperature and heat flux is proportional to the prescribed baseline error, which was to be expected. The 
results also demonstrate that the heat flux error bounds remain constant with time whereas the temperature error 
bounds decrease with time. Moreover, while the error in temperature stays well within the bounds of the perturbation 
(415%), the error in heat flux follows closely to the bounds of the perturbation. 


Sensitivity of Surface Solution to Error in Thermocouple Data 

A numerical study was also conducted to assess the sensitivity of the surface solution to the error in thermocouple 
data. For the sensitivity study, the embedded thermocouple is assumed to be at depth 31/4. 

To simulate error in thermocouple data, the stated methodology is solved with inaccurate temperature data located 
at accurate depths. The temperature data are perturbed by a factor of 6. For 6 = 0.01, 0.02, 0.03, 0.04, 0.05, 0.10, and 
0.15, the DHCP is solved with analytical solutions of T(x,,t) and T(x2,t) + 6T(x2,t) and depth inputs of x, = 0 
and x, = 3L/4. The solution yields qpycp, x, defined by Eq. (10). Boundary values T (x2, ¢) + 6T (x2, t) and qpucp, x, 
are then used to solve the IHCP to estimate the surface histories, T, and q,. 

Surface estimations are compared against the error in temperature and heat flux, e; and e,, respectively, obtained 
from the embedded thermocouple depth study. The known e7 and eg from depth x,2, shown in Figure 5 and Figure 
6, are used as a baseline in which case the perturbation, 6, is equal to 0.00. Since the model problem is a linear problem, 
the difference in error in the estimated temperature and heat flux is expected to be proportional to the prescribed error 
in the simulated baseline data. 

The Constant Case temperature and heat flux estimates are shown in Figure 9(a) and Figure 9(b) for At = 0.1 s. 
Additionally, the Linear Case temperature and heat flux estimates are shown in Figure 10(a) and Figure 10(b) for 
At = 0.1 s. Results demonstrate that for both the Constant and Linear Case, as 6 increases (decreases) towards +0.15 
(—0.15), the error in surface temperature and heat flux increase (decrease). Hence, the difference in error in the 
estimated temperature and heat flux is proportional to the prescribed baseline error, which was to be expected. The 
results also demonstrate that the temperature error bounds remain constant with time, whereas the heat flux error 
increases with time. Moreover, while the error in temperature follows closely to the bounds of the perturbation 
(+15%), the error in heat flux is well outside the bounds of the perturbation. 
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(b) Heat flux error at x = L 
Figure 7. Sensitivity of the surface estimations to the error in the specified location of the embedded 


thermocouple for the Constant Case using a time-step of At = 0.1 s. 
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Perturbation of embedded 
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(b) Heat flux error at x = L 
Figure 8. Sensitivity of the surface estimations to the error in the specified location of the embedded 


thermocouple for the Linear Case using a time-step of At = 0.1s. 
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Perturbation of measured 
temperature data, 6 
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Figure 9. Sensitivity of the surface estimations to the error in thermocouple data for the Constant Case using a 
time-step of At = 0.1 s. 
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temperature data, 6 
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(b) Heat flux error at x = L 


Figure 10. Sensitivity of the surface estimations to the error in thermocouple data for the Linear Case using a 
time-step of At = 0.1s. 
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IV. Concluding Remarks 

The purpose of this paper is to expand upon the 2016 work of Pizzo, et. al,! by characterizing the selected method 
for inverse heat conduction. Due to the ill-posedness of inverse problems, small errors in the data yield large errors in 
the solution. When solving inverse problems, accuracy in the estimated surface temperature and heat flux is dependent 
on several factors. Factors include the embedded thermocouple depth through the thickness of a material, the 
sensitivity of the surface solution to the error in the specified location of the embedded thermocouples, and the 
sensitivity to the error in thermocouple data. The effect of these factors on solution accuracy is presented in this paper. 
Each numerical study is performed using a model problem. The model assumes that one thermocouple is embedded 
through the thickness of a material and a second thermocouple is installed on the back face. Numerical data are 
compared against an exact analytical solution. 

Preliminary results of the study to determine if there is an optimal depth to embed one thermocouple through the 
thickness of a material demonstrate that the surface estimations for both temperature and heat flux are consistently the 
most accurate if one thermocouple is embedded between depths x19 and x12 (5L/8 and 3L/4) when a second 
thermocouple is installed on the back face. This finding is based entirely off of the model problem, and is independent 
of time-step. Further assessments should be conducted to conclusively make this assertion, by perhaps considering a 
radiative boundary condition for gq, with temperature-dependent thermophysical properties. This type of study would 
better emulate the testing conditions of carbon/carbon materials. Test data can then be used to validate the numerical 
findings. To assess whether the optimal depth is dependent on either the physical or numerical aspects of the problem, 
other spatial divisions should also be considered, e.g. by dividing the domain from x = 0 to x = L into perhaps 20 or 
100 evenly spaced intervals. 

When assessing the sensitivity of the surface solution to (a) the error in the specified location of the embedded 
thermocouple and to (b) the error in thermocouple data, the difference in error in the estimated temperature and heat 
flux is proportional to the prescribed baseline error, which was to be expected. The results also demontrate that when 
assessing the error in the specified location of the embedded thermocouple, the heat flux error bounds remain constant 
with time whereas the temperature error bounds decrease with time and when assessing the error in thermocouple 
data, the temperature error bounds remain constant with time, whereas the heat flux error increases with time. This 
disparity between sensitivity assessments is also to be expected. 

When perturbing the embedded depth, the thermocouple temperature remains unchanged from the baseline 
problem. In this case, the estimated surface temperature and heat flux are merely offset by a factor of the perturbation 
and no additional heat is added to, or removed from, the system. Hence, the results of depth perturbation follows 
closely to the baseline results of temperature [error decreases with time as indicated in Figure 5(a) and Figure 6(a)] 
and heat flux [error remains constant with time as indicated in Figure 5(b) and Figure 6(b)]. On the other hand, when 
perturbing the thermocouple data, additional heat is added to (+6) or removed from (—4) the baseline problem. With 
the addition (removal) of heat, the surface estimations for both temperature and heat flux will increase (decrease) with 
time as compared to the baseline problem. Hence, rather than the error in temperature decreasing with time as is the 
case for the baseline problem, the additional heat added/removed from the system causes the error to remain constant. 
Similarly, rather than the error in heat flux remaining constant with time as is the case for the baseline problem, the 
additional heat added/removed from the system causes the error to increase. This result demonstrates that any error in 
embedded depth will translate to a proportional error of the same order of magnitude, and any error in thermocouple 
data will translate to a proportional error of a higher order of magnitude. Therefore, when a user instruments a material 
with two thermocouples and post-processes the measured data to estimate the surface temperature and heat flux, it is 
of utmost importance that the measured depths and temperature data are accurate. 
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